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A zigzag graphene nanodisk can be interpreted as a quantum dot with an internal degree of 
freedom. It is well described by the infinite-range Heisenberg model. We have investigated its ther- 
modynamical properties. There exists a quasi-phase transition between the quasi-ferromagnet and 
quasi-paramagnet states, as signaled by a sharp peak in the specific heat and in the susceptability. 
We have also analyzed how thermodynamical properties are affected when two leads are attached 
to the nanodisk. It is shown that lead effects are described by the many-spin Kondo Hamiltonian. 
There appears a new peak in the specific heat, and the multiplicity of the ground state becomes just 
one half of the system without leads. Another lead effect is to enhance the ferromagnetic order. Be- 
ing a ferromagnet, a nanodisk can be used as a spin filter. Furthermore, since the relaxation time is 
finite, it is possible to control the spin of the nanodisk by an external spin current. We then propose 
a rich variety of spintronic devices made of nanodisks and leads, such as spin memory, spin amplifier, 
spin valve, spin-field-effect transistor, spin diode and spin logic gates such as spin-XNOR gate and 
spin-XOR gate. Graphene nanodisks could well be basic components of future nanoelectronic and 
spintronic devices. 



I. INTRODUCTION 

Graphene nanostructure^i^i^ has attracted much at- 
tention for its potential for future application in nano- 
electronics and spintronics. In particular, in graphene 
nanoribbons^i^i^i^^^'^'^^ the low-energy bands are almost 
flat at the Fermi level due to the edge states. Such a 
peculiar band structure has motivated many researchers 
to investigate their electronic and magnetic properties. 

Another basic element of graphene derivatives is a 
graphene nanodisk^^ ^^^i^^i^^i^^i^^i^^ . It is a nanometer- 
scale disk-like material which has a closed edge. There 
are many type of nanodisks, where typical examples are 
displayed in FigHJ Among them, the trigonal zigzag nan- 
odisk is prominent in its electronic and magnetic prop- 
erties because there exist A/'-fold degenerated half-filled 
zero-energy states when its size is N . Furthermore, spins 
make a strong ferromagnetic order due to the exchange 
interaction as large as the Coulomb interaction^^, and 
hence the relaxation time is finite but quite large even if 
the size N is very small. We have called such a system 
quasi-ferromagnet ^ ^ . 

In this paper we make a further study of quasi- 
ferromagnet by exploring thermodynamical properties of 
the nanodisk-spin system. The system is well described 
by the infinite-range Heisenberg model. A nanodisk is 
interpreted as a quantum dot with an internal degree of 
freedom. It is exactly solvable. Constructing the par- 
tition function, we calculate the specific heat, the en- 
tropy, the magnetization and the susceptibility. We find 
a sharp peak in the specific heat and in the susceptibility, 
which is interpreted as a quasi-phase transition between 
the quasi-ferromagnet and quasi-paramagnet states. 

We then investigate a nanodisk-lead system, where 
the lead is made of a zigzag graphene nanoribbon or an 
ordinary metallic wire. (We refer to it as a graphene 
lead or a metallic lead for brevity.) We perform the 
Schrieffer- Wolff transformation to derive the many-spin 
Kondo Hamiltonian describing the lead effects. Con- 
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FIG. 1: Basic configurations of typical graphene nanodisks. 
(a) Benzene, (b) Trigonal zigzag nanodisk (phenalene). (c) 
Trigonal armchair nanodisk (triphenylene). (d) Hexagonal 
zigzag nanodisk (coronene). (e) Hexagonal armchair nanodisk 
(hexa benzocoronene)^^. 



structing the partition function, we analyze thermody- 
namical properties. There appears a new peak in the 
specific heat but not in the susceptibility for small size 
nanodisks, < 10. Near the peak the internal energy is 
found to decrease. We show in the instance of the metal- 
lic lead that the band width of free electrons in the lead 
becomes narrower due to the Kondo coupling. We inter- 
pret this phenomenon to mean that some free electrons 
in the lead are consumed to make spin coupling with elec- 
trons in the nanodisk. Furthermore, the multiplicity of 
the ground state becomes just one half of that in the sys- 
tem without leads. They are indications of Kondo effects 
due to the Kondo interaction between electrons in the 
lead and the nanodisk. 

With respect to the ferromagnetic order, we find that 
the lead effect is to enhance the order. This is an impor- 
tant property to fabricate spintronic circuits by connect- 
ing leads to nanodisks in nanodevices. Being a ferromag- 
net, a nanodisk can be used as a spin filter. Furthermore, 
since the relaxation time is finite, it is possible to con- 
trol the spin of the nanodisk by an external spin current. 
We propose a rich variety of spintronic devices made of 
nanodisks and leads, such as spin memory, spin amplifier, 
spin valve, spin-field-effect transistor, spin diode and spin 
logic gates such as spin-XNOR gate and spin-XOR gate. 
Graphene nanodisks could well be basic components of 
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FIG. 2: (a) Geometric configuration of trigonal zigzag nan- 
odisks. It is convenient to introduce the size parameter N in 
this way. The 0-trigonal nanodisk consists of a single Ben- 
zene, and so on. The number of carbon atoms are related as 
Nc = N'^ + 6N + 6. (b) Density of states of the A^-trigonal 
nanodisk for N = 0, 1, 2, • • • ,7. The horizontal axis is the 
size N and the vertical axis is the energy e{N) in units of 
t — 2.7eV. Dots on colored bar indicate the degeneracy of 
energy levels. 

future nanoelectronic and spintronic devices. 

This paper is organized as follows. In SecUll we sum- 
marize the basic notion of trigonal zigzag nanodisks. The 
low-energy physics is described by electrons in the A^- 
fold degenerated zero-energy sector, which form a quasi- 
ferromagnet due to large exchange interactions. We also 
analyze thermodynamical properties of the nanodisk-spin 
system. In Sec lIIIt we derive the many-spin Kondo 
Hamiltonian by the Schrieffer- Wolff transformation in the 
nanodisk-spin system coupled with graphene leads and 
also metallic leads. The partition function is calculated 
in the two steps. First, we perform a functional inte- 
gration over the lead-electron degree of freedom. Sec- 
ond, we sum up over the nanodisk-spin degree of free- 
dom. We then analyze thermodynamical properties of 
the nanodisk-lead system. In Sec HVl we propose some 
spintronics devices made of nanodisks and leads. 

II. NANODISK QUASI-FERROMAGNETS 
A. Zero-Energy Sector 

We calculate the energy spectra of graphene deriva- 
tives based on the nearest-neighbor tight-binding model, 
which has been successfully applied to the studies of car- 
bon nanotubes and nanoribbons. The Hamiltonian is 
defined by 

H = Y,e4c,^Y.^,jc\cj, (2.1) 

where £i is the site energy, tij is the transfer energy, and 
ct is the creation operator of the tt electron at the site 
i. The summation is taken over all nearest neighbor- 
ing sites j). Owing to their homogeneous geometrical 
configuration, we may take constant values for these en- 
ergies, £i = Sy and tij = t ^ 2.70eV. Then, the diagonal 
term in (|2.ip yields just a constant, epNc, where Nq is 



the number of carbon atoms in the system. The Hamil- 
tonian (|2.ip yields the Dirac electrons for graphene^-i^i^. 
There exists one electron per one carbon and the band- 
filling factor is 1/2. It is customary to choose the zero- 
energy level of the tight-binding Hamiltonian (|2.ip at this 
point so that the energy spectrum is symmetric between 
the positive and negative energy states. Therefore, the 
system is metallic provided that there exists zero-energy 
states in the spectrum. A comment is in order. It is 
understood that carbon atoms at edges are terminated 
by hydrogen atoms. We carry out the calculation^'^^ to- 
gether with this condition. 

It is straightforward to derive the energy spectrum Ei 
together with its degeneracy gi for each nanodisk by di- 
agonalizing the Hamiltonian (|2.ip . The density of state 
is given by 

D{e) = Y^g,5{e-E,). (2.2) 

i=l 

We have found that the emergence of zero-energy states is 
surprisingly rare. Among typical nanodisks, only trigonal 
zigzag nanodisks have degenerate zero-energy states and 
show metallic ferromagnetism. As an example we display 
the density of state (|2.2p of trigonal zigzag nanodisks in 
FiglS We have classified them by the size parameter N 
as defined in Figl2l(a), in terms of which the number of 
carbons is given by Nq = N'^ + QN + 6. 

The size-TV nanodisk has A^-fold degenerated zero- 
energy states, where the gap energy is as large as a 
few eV. Hence it is a good approximation to investigate 
the electron-electron interaction physics only in the zero- 
energy sector, by projecting the system to the subspace 
made of those zero-energy states. The zero-energy sector 
consists of N orthonormal states \fa)j <^ = 1, 2, • • • , A^, 
together with the SU(A^) symmetry. We can expand the 
wave function of the state \fa) as 

Ux) = Y,^?^^i^)^ (2.3) 

i 

where (pi{x) is the Wannier function localized at the site 
i. The amplitude a;f is calculable by diagonalizing the 
Hamiltonian (|2.ip . All of them are found to be real. It 
is intriguing that one of the wave functions is entirely lo- 
calized on edge sites for the nanodisk with TV =odd, as in 
FiglSl where the solid (open) circles denote the amplitude 
uji are positive (negative) . The amplitude is proportional 
to the radius of circle. Such a wave function does not ex- 
ist for the nanodisk with N =even. 



B. Coulomb Interactions 

We include the Coulomb interaction between electrons 
in the zero-energy sector^-. It is straightforward to 
rewrite the Coulomb Hamiltonian as Hj^ = Hs -\- Hu 
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FIG. 3: The zero-energy states of the trigonal nanodisk with 
size N = 5. There are 5 degenerate states. The soUd (open) 
circles denote the amplitude Ui are positive (negative). The 
amplitude is proportional to the radius of circle. Electrons 
are entirely localized on edges in one of the states. 



with 



(2.4a) 



a>f3 



= [Uaf3 - ^-^«/3 j n (a) n (/?) + ^ Uaan {a) , 

(2.4b) 



a>f3 



where Ua/3 and Ja/3 are the Coulomb energy and the ex- 
change energy between electrons in the states \fa) and 
I//3). Here, n {a) is the number operator and S{a) is the 
spin operator, 

n{a) = ^^^dl^{a)da{a), (2.5a) 

a 

S{a) = i^4(a)r,,.d,.(a), (2.5b) 



with da (a) the annihilation operator of electron with spin 
cr =t, I in the state \fa)', 't is the Pauli matrix. 

The remarkable feature is that there exists a large over- 
lap between the wave functions /^(ic) and //3(cc), 7^ 
since the state \fa) is an ensemble of sites as in (|2.3p 
and identical sites are included in \foi) and I//3). Conse- 
quently, the dominant contributions come from the on- 
site Coulomb terms not only for the Coulomb energy 
but also for the exchange energy. Indeed, it follows that 
Ua(3 = Ja(3 in the on-site approximation. We thus obtain 



(2.6) 



where 



U = jd^xd^y^;{x)^,{x)V{x-y)^;{y)^,{y), (2.7) 

with the Coulomb potential V{x — y). The Coulomb en- 
ergy U is of the order of leV because the lattice spacing of 
the carbon atoms is ~ lA in graphene. Coulomb block- 
ade peaks appear (i = and /i = Sq, + Ucx(3^ where 
new channels open^^. We can determine experimentally 
the energy £0, and Ua,(3 by identifying Coulomb blockade 
peaks. 



C. SU(A/^) Approximation 

Since the exchange energy 7^/3 is as large as the 
Coulomb energy /7q,/3, the spin stiffness Ja(3 is quite large. 
Furthermore, we have checked^ ^ numerically that all 7^/3 
are of the same order of magnitude for any pair of a and 

implying that the SU(A^) symmetry is broken but not 
so strongly in the Hamiltonian (|2.4ap . It is a good ap- 
proximation to start with the exact SU(A/') symmetry. 
Then, the zero-energy sector is described by the Hamil- 
tonian i^D = + , with 



jY,S{a)-S{(3), (2.8a) 
^ ^n(a)n(/?) + /7^n(a), (2.8b) 



where J ^ U. The term Hs is known as the infinite-range 
Heisenberg model. We rewrite them as 



JS 

U 
2" 



tot 

J 

4 



3 J 

K'ot + 1), 



(2.9a) 
(2.9b) 



where ^Stot = (^) total spin, and ntot = 

n (a) is the total electron number. 
The ground states of nanodisks are half filled. We re- 
strict the Hilbert space to the half-filling sector. 



n {a) = {a) + {a) = 1. 



(2.10) 



The Hamiltonian (|2.9ap is reduced to the Heisenberg 
model. 



Hs 



-J Stot • Sf, 



(2.11) 



where we have neglected an irrelevant constant term, 
(3/4) J7V. This is exactly diagonalizable, i^sl^) = 
Esl"^), with 



-J5(5+ 1), 



(2.12) 



where s takes values from N/2 down to 1/2 or 0, depend- 
ing on whether N is odd or even. 



N N 



1, 



N 



2,- 



2 ' 2 '2 
The Hilbert space is diagonalized. 



s > 0. 



^1/2 



(2.13) 



(2.14) 



where Hs denotes the {2s + 1) dimensional Hilbert space 
associate with an irreducible representation of SU(2). 
The multiplicities {s) satisfies the recursion relation 
coming from the spin synthesizing rule. 



9n {s) qn- 



+ gN-i I s + 2 



(2.15) 
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We solve this as 



(a) 1.2 



N 



9N[—-q 



Cq —nCq-I. 



(2.16) 



IS 



The total degeneracy of the energy level E, 
{2s^l)gN{s). 

At half filling, the eigenstate of the Hamiltonian Hj^ is 
labeled as |^) = Intot, where s is the total spin and 

m is its component. We refer to the total spin 5'tot of 
a nanodisk as the nanodisk spin. 



D. Thermodynamical Properties 

We have a complete set of the eigenenergies together 
with their degeneracies. The partition function of the 
nanodisk with size N is exactly calculable, 



Zs = Y,{2s^l)gN{s)e- 



N/2 

J2{N-2q+l) {nC, 

9=0 



N Cq-l) 



X exp 



(2.17) 



According to the standard procedure we can evaluate the 
specific heat C{T), the entropy S{T), the magnetization 

(-^tot) and the susceptibility X = {{^tot) - i^totf) 
from this partition function, where 



(2.18) 



is the ground-state value of the total spin. The entropy 
is given by 




S (0) = kB log(iV + 1) 



(2.19) 



at zero temperature. We display them in FigHlfor size 
A/' = 1,2,22,- • -210^ 

There appear singularities in thermodynamical quan- 
tities as ^ oo, which represent a phase transition at 
Tc between the ferromagnet and paramagnet states. 



JN 
2k^' 



(2.20) 



For finite there are steep changes around Tc, though 
they are not singularities. It is not a phase transi- 
tion. However, it would be reasonable to call it a 
quasi-phase transition between the quasi-ferromagnet 
and quasi-paramagnet states. Such a quasi-phase transi- 
tion is manifest even in finite systems with N = 100 ~ 
1000. 

The specific heat and the magnetization take nonzero- 
values for T > Tc [Figl4Ka),(c)], which is zero in the limit 




FIG. 4: Thermodynamical properties of the nanodisk-spin 
system, (a) The specific heat C in unit of /cb^V. (b) The 
entropy S in unit of A:BiVlog2. (c) The magnetization (S'tot) 
in unit of Sg. (d) The susceptibility x in unit of Sg. The 
size IS iV = 1, 2, 2^ • • • 2^°. The horizontal axis stands for the 
temperature T in unit of JN/kB- The arrow represents the 
phase transition point Tc in the limit N ^ oo. 



nanodisk-lead system 




FIG. 5: The nanodisk-lead system. The nanodisk with N = 7 
is connected to the right and left leads by tunneling coupling 
tR and tL- 



^ oo. The entropy for T > Tc is lower than that 
of the paramagnet [FigH^b)]. These results indicate the 
existence of some correlations in the quasi-paramagnet 
state. The maximum value of the susceptibility increases 
linearly as N becomes large. It is an indicator of the 
quasi-phase transition. 



III. MANY-SPIN KONDO EFFECTS 



We proceed to investigate how thermodynamical prop- 
erties of the nanodisk is affected by the attachment of 
the leads. The nanodisk is no longer in the half-filled 
state, when charges transfer between the nanodisk and 
the leads. However, the nanodisk remains to be half 
filled, when a charge transfers from the lead to the nan- 
odisk and then transfers back from the nanodisk to the 
lead. Such a process is the second order effect in the 
tunneling coupling constant i. We now show that it is 
described by the many-spin Kondo Hamiltonian. 
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A. Nanodisk-Lead System 



B. Many-Spin Kondo Hamiltonian 



We analyze a system made of a nanodisk connected to 
two leads [FiglS]. The model Hamiltonian is given by 



Here, i^L is the lead Hamiltonian, 



(3.1) 



(3.2) 



ka 



describing a noninteracting electron gas in the leads with 
the dispersion relation e (k), where c^^ is the annihilation 
operator of electron with the wave number k and the 
spin a in the left (x =L) or right (x =R) lead. On the 
other hand, H^l and H^r are the transfer Hamiltonians 
between the left (L) and right (R) leads and the nanodisk, 
respectively. 



ffxL =tL E + 4(«)cL], (3-3a) 

ka a 

Htr =tR Y.T.^^'klMa) + di{a)cfj, (3.3b) 

ka a 

with the tunneling coupling constant. We have as- 
sumed that the spin does not flip in the tunneling pro- 
cess. 

The nanodisk-lead system looks similar to that of the 
A^-dot system. However, there exists a crucial differ- 
ence. On one hand, in the ordinary A/'-dot system, an 
electron hops from one dot to another dot. On the other 
hand, in our nanodisk system, the index a of the Hamil- 
tonian runs over the A'-fold degenerate states and not 
over the sites. According to the Hamiltonian (|3.3p . an 
electron does not hop from one state to another state. 
Hence, it is more appropriate to regard our nanodisk as 
a one-dot system with an internal degree of freedom. 

It is convenient to make the transformation 



^ka 
^ka 



1 



-^ka 
^R 

■^ka 



(3.4) 



with 



so that the right and left leads are combined into the 
"even" and "odd" leads. The lead Hamiltonian Hi, is 
invariant under above transformation, 

HL = ^e ik) (4cL + 4C) , (3.6) 

ka 

but the transfer Hamiltonian is considerably simplifled, 
= i E E (4lrf<.(«) + 4(«)cL) ■ (3-7) 

ka OL 



It looks as if the tunneling occurs only between the 
"even" lead and the nanodisk. Namely, we may neglect 
the "odd" lead term in the lead Hamiltonian (|3.6p . 



The total Hamiltonian is = i^s + i^u + + ^t- 
We analyze the Hamiltonian H = Hq ^ H^j by taking 
Ho = Hu + i^L as the unperturbed term and Ht as 
the perturbation term. We make a canonical transfor- 
mation known as the Schrieffer- Wolff transformation^, 
H ^ H = e^^i^e"*^, with G the generator satisfying 



Ht + -[G,Ho]=0. 



(3.8) 



We may solve this condition explicitly for the generator, 

i 



G 



-iE( 

Ni [^\ 



€d-£ (k) 
i 



_ rf ef 
^aa^ka^<^^ 



e{k')- 

ee{ 

Oika (3a' 



t 



e (k') - £d 



ni3o 



h.c. 



where cr =|t for a =t|, and 

Ja(3 



Ua(3 ■ 



(3.9) 



(3.10) 



The leading term is the second order term, and given by 



H^ 



(2) 



eff 



27V 



a(3kk' a 



N ^ 



ed-e{k)^U'^^ e{k')-ed ) 
1 1 



-^k' a 



a(3kk' aa' 



X c 



Taa'cl,^, ■S{a). 



(3.11) 



The dominant contribution comes from the Fermi sur- 
face, £(k) = sf- We now assume the SU(A^) symmetry 



(3.5) U^^Q = and the symmetric condition sp = Sd- 



u' 



with 

respect to the Fermi level. Then, the second order term 
becomes the many-spin Kondo Hamiltonian, 



Hy. = i^eff = ^ Ck^r aa'Cl^a' ' 

kk' aa' 

with the Kondo coupling constant 
1 1 



(3.12) 



Jk = 2t' 



£d - ep + U' ep - £d 



8*2 
U' 



777- (3-13) 



The difference between the above many-spin Kondo 
Hamiltonian and the ordinary Kondo Hamiltonian is 
whether the local spin is given by the summation over 
many spins ^Stot or a single spin S. Note that S'^q^ is a 
dynamical variable but S'^ is not, S'^ =3/4. 
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C. Kondo-Heisenberg System 



Performing the integration we find 



We have derived the Kondo Hamiltonian from the 
Coulomb and transfer terms by way of the Schrieffer- 
Wolff transformation. The resultant system is the 
Kondo-Heisenberg model, which comprises of the Heisen- 
berg Hamiltonian, the lead-electron Hamiltonian and the 
Kondo Hamiltonian, 



Zk = Det[M] =exp [-(3Fk].. 



(3.21) 



(3.14) 



where we have ignored Hjj since it is just a constant at 
half filling. We note that the order of the Heisenberg 
Hamiltonian (~ U) is much larger than the order of the 
Kondo Hamiltonian (~ 4P//7) because U ^t. 

Our goal is the analysis of the partition function of the 
coupled system (|3.14p . We define the spinor ^jj — (cp cp^ 
The partition function in the Matsubara form is given by 



Z = Tr^ / Vi)Vi)^ exp 



(3.15) 

with TYeff the Hamiltonian density, where the lead elec- 
tron's degree of freedom is integrated out by a functional 
integral, and the nanodisk spin is summed up. Since the 
Heisenberg term does not contain lead electrons, we can 
separate it as 

Z = TTs[e^v{-pHs)Z^], (3.16) 
where Zk = J VipVi/j^ exp [— 6'k], with 

Sk = j ii^^dri^ ^Hl^ Hk) . (3.17) 

First, we evaluate the functional integration Zk in Sub- 
section [IIIdI We obtain the effective spin Hamiltonian, 
where the only active degree of freedom is the nanodisk 
spin. Then, we sum up over the nanodisk spin to ob- 
tain the total partition function Z. We do this for the 
graphene lead and for the metallic lead in Subsections 
nil El and IIII F[ respectively. 



D. Functional Integration 

Because an electron in the lead is constrained within a 
very narrow region, it is a good approximation to neglect 
momentum scatterings. 



kcTcr' 

The action (|3.17p is summarized as 

/ ^^V^(fc)M(fc)^(fc), 

k 

with 

M (k) = -[ico - e {k)] + Jkt • S'tot- 



(3.18) 



(3.19) 



(3.20) 



where 



i^K = - ^ ^ In [cosh/3JK l^totl + cosh/3e {k)] . (3.22) 



is the Helmholtz free energy Fk. This formula is reduced 
to the well-known one for free electrons with the disper- 
sion relation e{k) for Jk = 0. We evaluate the momen- 
tum integral for the graphene lead and for the metallic 
lead, separately, in the succeeding two subsections. 



E. Zigzag Graphene Nanoribbon Leads 

We consider the system where the leads are made of 
zigzag graphene nanoribbons. Owing to the fiat band at 
the zero energy, e (k) = 0, the result of the functional 
integration (|3.22p is quite simple. 



Fk = -^Incosh^Jx |5'tot| 



(3.23) 



The effective Hamiltonian for the nanodisk spin is Hs + 
Fk. The lead effect is to make the effective spin stiffness 
larger and the ferromagnet more rigid. 
The partition function (|3.16p is 



Z = Tts 



exp [/^J^Stot] cosh ^ Jk \S- 



(3.24) 



as implies that the eigenstates of the total Hamiltonian 
are given by JS'^^^ ± ^ Jk \ 'S'tot | • Namely, the ground 
states are split into two by the existence of the lead- 
electron spin. Accordingly, the ground state multiplicity 
is reduced to (A^ + 1) /2, which is just one half of that of 
the nanodisk without leads. 

The trace over the total spin is carried out in (|3.24p . 



N/2 

Z = Y,{N-2q^l) {NCq -N Cq- 

q=0 

'TV 



X exp 
X cosh 



'N \ [N 



. (3.25) 



In FiglSl we show the specific heat Cq{T)^ the entropy 
Sg{T), the magnetization ('S'^ot) susceptibility 
X for various size N connected with graphene leads. 

We compare thermodynamical properties of the nan- 
odisk without leads (FigS]) and with leads (FigjG]). The 
significant feature is the appearance of a new peak in the 
specific heat at Tk = (Jk/2J)Tc, though it disappears 
for large N. We examine the internal energy Eq{T), 
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FIG. 6: Thermo dynamical properties of the nanodisk-spin 
system with graphene leads, (a) The specific heat C in unit 
of ksN. (b) The entropy S in unit of (A:BiVlog2). (c) The 
magnetization (/Stot) in unit of Sg. (d) The susceptibility 
X in unit of Sg. The size is = 1, 2, 2^ • • • 2^°. We have 
set Jk/J = 0.2. The horizontal axis stands for the temper- 
ature T in unit of JN/k^- The arrows represents the points 
corresponding to Tc and Tk. 



which is found to decrease around Tk (FigE)). Near zero 
temperature it reads 



Eg{T) 



-J St 



Jk 



Sg + J^Sge 



(3.26) 



The first term represents the energy stabilization due to 
the ferromagnetic order present in the nanodisk system 
without leads, while the second term represents the one 
due to the Kondo coupling Jk between spins in the nan- 
odisk and in the leads. Furthermore, it follows that the 
entropy is reduced at zero temperature as 

5G(0)-^(0) = -felog2 (3.27) 

with (|2.19p . as implies that the ground state multiplicity 



(a) without leads 

0.2 0.4 ^ 0.6 



(b) with graphene leads 

0.2 



-0.2 




-0.8 



FIG. 7: The internal energy E/NkB for the system (a) 
without leads and (b) with graphene leads, for size N — 
1, 2, 4, • • • 1024. The horizontal axis stands for the tempera- 
ture, (a) The energy descreases except for = 1 as the tem- 
perature decreases, which represents the ferromagnetic order, 
(b) There exists an additional energy decrease around Tk, 
which is prominent for = 1, attributed to the Kondo effect. 
The arrows represents the points corresponding to Tc and Tk. 



at the zero temperature is just one half of that of the sys- 
tem without leads, in accord with the observation made 
below (|3.24p . These features indicate the occurrence of 
the Kondo effect due to the coupling between the spins 
in the nanodisk and the leads. 



F. Metallic Leads 

Next we consider the system comprised of metallic 
leads with a constant energy density. 
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\e\<D 
\e\>D 



(3.28) 



We change the momentum integration into the energy 
integration in (|3.22p . 

/D 
de In [cosh/?JK l^totl + cos\i(3e\ . (3.29) 
-D 

The integration is carried out as 



PDJ^\Stot\-2D\og2 

/3(D-JK|Stot|)^ 

^-/3(D+JK|Stot|)^ 



-Li3(-e-«°+*l«-l')} 
where Li2 [x] is the dilogarithm function—, 

[A = Y.^' 



(3.30) 



(3.31) 



It is easy to see that the free energy is reduced to that 
of the nanodisk-spin system with graphene leads in the 
limit T> ^ with p 1/2D. 

In FiglH we show the specific heat Cm(T), the entropy 
Sm(T)^ the magnetization (^S^ot) susceptibility 
X for various size N connected with metallic leads. 

The behaviors of the magnetization and the suscepti- 
bility are the same as those in the case of the nanodisk- 
spin system with graphene leads. On the other hand, 
overall behaviors are the same with respect to the specific 
heat and the entropy. In particular, the same relation as 
(|3.27p holds for the entropy. 



5M(0)-5(0) = -feB log 2. 



(3.32) 



There are some new features in low temperature regime. 
Using the asymptotic behavior^ 



lim Li2 [- 



lim Lio I 



TT 



1. 2I 

2 X 



0, 



(3.33) 



we obtain the free energy, the entropy, the specific heat 
and the internal energy up to the terms in the order of 
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FIG. 8: Thermo dynamical properties of the nanodisk-spin 
system with metaUic leads, (a) The specific heat C in unit 
of kBN. (b) The entropy S in unit of kBN\og2. (c) The 
magnetization {S^^^) in unit of S'^. (d) The susceptibility x 



in unit of Sg. The size is = 1, 2, 2^ 



We have set 



Ji^/J — 0.2 and D — 2/cbTc. The horizontal axis stands for 
the temperature T in unit of JN/kB. The arrows represents 
the points corresponding to Tc and Tr. 



as follows, 



(T) 
Cm{T) 

and 



kB log 

^2 



tot I 



2Dlog2 + -— 



TT 

Op 



IV. SPINTRONIC DEVICES AND 
APPLICATIONS 

We propose some applications of graphene nanodisk- 
lead systems to spintronic devices^^. The nanodisk-spin 
system is a quasi-ferromagnet, which is an interpolating 
system between a single spin and a ferromagnet. It is 
easy to control a single spin by a tiny current but it 
does not hold the spin direction for a long time. On the 
other hand, a ferromagnet is very stable, but it is hard to 
control the spin direction by a tiny current. A nanodisk 
quasi-ferromagnet has an intermediate nature: It can be 
controlled by a relatively tiny current and yet holds the 
spin direction for quite a long time. Indeed, its life-time 

Tferro IS givCU hy^ 



Tferro OC exp 



2kT 



(4.1) 



which is quite long compared to the size. The important 
point is that the size is of the order of nanometer, and it 
is suitable as a nanodevice. 

The coupling of the nanodisk spin and the in- 
jected electron spin ^i/jWiIj is described by the Landau- 
Lifschitz-Gelbert equation^. 



TT 

~3 



pklT, 



- Jk l-^totl) 

(3.34) 
(3.35) 

(3.36) 



dn 

— 7Beff X n 
at 



an X 



dn 



(4.2) 



where n = ^Stot / 1 '^tot | is the normalized nanodisk spin, 
7 is the gyromagnetic ratio, a is the Gilbert damping 
constant {a ^ 0.01), and ^eff is the effective magnetic 
field produced by the injected electron spin. 



U 



>eff ^ 



tot I 



(4.3) 



Em{T)=Eg (T) + A^(T). 



(3.37) 



It is proportional to the injected current P^. A spin 
polarized current rotates the nanodisk spin to the same 
direction as the current with the relaxation time 



The specific heat Cm {T) is identical to the specific heat 
of free electrons in the metallic lead. The internal energy 
Em (T) consists of two terms: Eq (T) is identical to the 
energy (|3.26p for the nanodisk with graphene leads, and 



Tfilter = 



oc A^. 



We use these properties to design spintronic devices. 



(4.4) 



AEiT) --^{D- JKSgf + y p {ksTf (3.38) 

is the energy for the metallic lead. Here, the first term 
shows that the band width of free electrons in the lead 
(|3.28p becomes narrower due to the Kondo coupling. We 
may interpret that n free electrons in the lead with 

n = pJ^Sg (3.39) 

are consumed to make spin-coupling with electrons in the 
nanodisk. The second term is the thermal energy of free 
electrons in the metallic lead. 



A. Basic Components of Spintronic Devices 

Spin filter: We consider a lead-nanodisk-lead system 
[Figll], where an electron makes a tunnelling from the 
left lead to the nanodisk and then to the right lead. Lead 
electrons with the same spin direction as the nanodisk 
spin can pass through the nanodisk freely. However, 
those with the opposite direction feel a large Coulomb 
barrier and are blocked (Pauli blockade)^. As a result, 
when we apply a spin-unpolarized current to the nan- 
odisk, the outgoing current is spin polarized to the di- 
rection of the nanodisk spin. Consequently, this system 
acts as a spin filter. 



Spin memory: A nanodisk can be used as a spin mem- 
ory, where the spin direction is the information. We can 
read-out the information by applying a spin-unpolarized 
current because the outgoing current from a nanodisk is 
spin-polarized to the direction of the nanodisk spin. Fur- 
thermore, the direction of the nanodisk spin itself can be 
controlled by applying a spin-polarized current into the 
nanodisk. 

Spin amplifier: A nanodisk can be used as a spin 
amplifier. We take the incoming current to be par- 
tially polarized, whose average direction is assumed to 
be up, lY > I'l" > 0. On the other hand, the di- 
rection of the nanodisk spin is arbitrary. Since spins 
in the nanodisk feel an effective magnetic field propor- 
tional to — /|", they are forced to align with that of 
the partially-polarized-spin current after making damped 
precession. After enough time (r Tfiiter), all spins in 
the nanodisk take the up direction and hence the outgo- 
ing current is the perfectly up-polarized one, I^^^ = P^, 



spin valve 



rout 



0. Consequently, the small difference IV 



IS 



amplified to the large current . The amplification ratio 
is given by P^/ {P^ — /|"), which can be very large. This 
effect is very important because the signal of spin will 
easily suffer from damping by disturbing noise in leads. 
By amplifying the signal we can make circuits which are 
strong against noises. 

Spin rotator: We can arrange a lead so that it has the 
Rashba-type interaction^^. 



(4.5) 



Spins precess while they pass through the lead. The spin- 
rotation angle is given^^ by ^ = 2Am*L/fi, where m* is 
the electron effective mass in the lead and L is the length 
of the lead. We can control by changing the coupling 
strength A externally by applying an electric field^^. In 
this way we can rotate the direction of spin current by 
any degree 0. We call such a lead as a spin rotator. 



B. Some Spintronic Devices 

Spin valve: A nanodisk can be used as a spin valve, 
inducing the giant magnetoresistance effec t^^i^^i^^ . We 
set up a system composed of two nanodisks sequentially 
connected with leads [Figl9]. We apply external magnetic 
field, and control the spin direction of the first nanodisk 
to be \0) — cos I It) + sin I II), and that of the second 
nanodisk to be |0) = |T)- We inject an unpolarized-spin 
current to the first nanodisk. The spin of the lead be- 
tween the two nanodisks is polarized into the direction 
of |6>). Subsequently the current is filtered to the up-spin 
one by the second nanodisk. The outgoing current from 

/ cos I . We can control 



the second nanodisk is 1^^^ 

the magnitude of the up-polarized current from to / by 
rotating the external magnetic field. The system act as 
a spin valve. 




(b) 



//// 



R 



FIG. 9: Illustration of spin valve, (a) The spin valve is made 
of two nanodisks with the same size, which are connected 
with leads, (b) Applying external magnetic field, we control 
the spin direction of the first nanodisk to be |^), and that of 
the second nanodisk to be |0) = |T)- The incomming current 
is unpolarized, but the outgoing current is polarized, /|^^ = 
0. 



/cosf, /f^^ 



(a) 



spin-field-effect transisitor 





11 f My M 
i \ 

FIG. 10: Illustration of spin-field-effect transistor, (a) It is 
made of two nanodisks with the same size, which are con- 
nected with a rotator, (b) We set the spin direction of the 
two nanodisks to be up by magnetic field. The incomming 
current is unpolarized, but the outgoing current is polarized 
and given by 1^"^^ = /cosf, /J''^ = 0. The up-spin current 
is rotated by the angle Q within the central lead ascting as a 
rotator. 



Spin- field- effect transistor: We again set up a system 
composed of two nanodisks sequentially connected with 
leads [FigHO]. We now apply the same external magnetic 
field to both these nanodisks, and fix their spin direction 
to be up, |0) = It). As an additional setting, we use a lead 
acting as a spin rotator with the spin-rotation angle 0. 
The outgoing current from the second nanodisk is I^^^ = 

/cos|. It is possible to tune the angle by applying 
an electric field. Hence we can control the magnitude of 
the up-polarized current. The system acts as a spin-field- 
effect transistor^^. 

Spin diode: We set up a system composed of two nan- 
odisks sequentially connected with leads, where two nan- 
odisks have different sizes [FigHT]. The left nanodisk is 
assumed to be larger than the right nanodisk. Then the 
relaxation time of the left nanodisk r^i^^^ is larger than 
that of the right nanodisk r^i^^^, '^mter > '^mter- Second, 
we apply the same magnetic field to the two nanodisks, 
but we take it so small that the nanodisk spin can be 
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spin diode 
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(b) 




!) spin rotator 
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FIG. 11: Illustration of spin diode, (a) It is made of two nan- 
odisks with different size, which are connected with a rotator, 
(b) We set initially the spin direction of the two nanodisks to 
be up by weak magnetic field. By controlling the bias voltage, 
the current flows from the left lead to the right lead, or in the 
opposite way. (c) We inject an unplolarized current from left 
(right), which is made up-polarized by the left (right) nan- 
odisk, and then rotated by the central lead. The spin of the 
right (left) nanodisk is rotated after the relaxation time r^^ter 
("^filter)- When the incoming current is an unpolarized pulse, 
the charge transported by the current is different whether it 
is injected from left or right. This acts as a spin diode. 



controlled by a polarized current. For definiteness we 
take the direction of the magnetic field to be up. Third, 
the central lead is taken to be a spin rotator with ^ ir. 
When no currents enter the nanodisk, the directions of 
the two nanodisk spins are identical due to the tiny ex- 
ternal magnetic field, which is up. This is the "off" state 
of the spin diode. When we inject an unpolarized cur- 
rent to this system, the outgoing current is initially very 
small, /°^* = I cos 6 c:^ for AO ^ tt. However, after the 
relaxation time Tfiiter, the outgoing current becomes large 
since the polarized current rotates the spin of the second 
nanodisk by the angle 6. It takes a longer time when the 
nanodisk size is larger. Now, let us inject an unpolarized 
pulse current either from left or from right. It follows 
that Q^^^ > Q^^^ since r^^^^^ > r^^^^, where Q^^^ 
and Q^^^ are the charges transported by the current 
injected to the right and left leads, respectively. Hence 
the system acts as a spin diode. 



C. Spin Logic Gates 

We can construct spin logic gates in which the spin 
direction takes logic values; truth (false) identified with 
up (down) spin, by controlling a spin current by another 
spin current according to the following setups. 

Spin inverter (spin NOT gate): We take a spin rotator 
with the rotation angle = ir. This rotator is used as a 
spin NOT gate, by regarding up spin as "true" and down 
spin as "false", because it interchanges up spin with down 
spin. 



(a) spin-XNOR gate 




(b) parallel external fields (spin current flows) 





(c) anti-parallel external fields (spin current is blocked) 




FIG. 12: (a) Illustration of spin-XNOR gate with unpolarized 
current coming from left. The unpolarized current is filtered 
by the left nanodisk, and only up-spin electrons go through 
the centeral lead, (b) When we apply vertical spin currents 
with parallel spin direction, the out going current exist, (c) 
When we apply vertical spin currents with antiparallel spin 
direction, the out going current does not exist. 



Spin XNOR: We may construct a spin XNOR gate. 
We set up a system composed of two nanodisks connected 
with seven leads, three horizontal leads and four vertical 
leads, as illustrated in Fig. [121 We control the nanodisk 
spins by vertically applied spin currents instead of ex- 
ternal magnetic field. We inject an unpolarized current 
from the left lead to the right lead. When the spin di- 
rections of the two vertical spin currents are parallel, the 
directions of the two nanodisk spins become parallel, and 
the horizontal current can pass through. However, when 
the spin directions of the two vertical spin currents are 
antiparallel, the horizontal current cannot pass through. 
This system act as a spin XNOR gate by regarding up 
spin as "true" and down spin as "false". 

Spin XOR: We may construct a spin XOR gate by 
changing the central lead to be a spin rotator with the 
rotation angle A6> = tt. The horizontal current passes 
through if the spin directions of the two vertical spin 
currents are antiparallel and cannot pass through if they 
are parallel, as illustrated in Fig. [131 
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(a) spin-NOR gate 
lULl 



UR 




|DL| |DR| 

(b) parallel external fields (spin current is blocked) 



(c) anti-parallel external fields (spin current flows) 




FIG. 13: (a) Illustration of spin-XOR gate with unpolarized 
current coming from left. The unpolarized current is filtered 
by the left nanodisk, and only up-spin electrons go through 
the centeral lead. The electron spin in the central lead ro- 
tates ^ = TT by the Rashba-type interaction, (b) When we 
apply vertical spin currents with parallel spin direction, the 
out going current does not exist, (c) When we apply vertical 
spin currents with antiparallel spin direction, the out going 
current exist. 



V. CONCLUSIONS 

The trigonal zigzag nanodisk has a remarkable prop- 
erty that it has A^-fold degenerate zero-energy states with 
the S\J{N) symmetry when its size is TV. The S\J{N) 
symmetry is broken but not so strongly by the Coulomb 
interactions. The system is well approximated by the 



infinite-range Heisenberg model, where the site index 
runs over these A^-fold degenerate states. We may regard 
it as a quasi-ferromagnet characterized by the exchange 
energy as large as the Coulomb energy. The relaxation 
time is finite but quite large even if the size is very small. 

In this paper we have investigated thermodynamical 
properties of a nanodisk. We have found the emergence of 
a quasi-phase transition between the quasi-ferromagnet 
and the quasi-paramagnet as a function of temperature 
even for samples with N ^ 100. The transition point Tc 
is signaled by a sharp peak in the specific heat and in the 
susceptibility. 

We have also examined how they are modified when 
the external leads are attached to the nanodisk. The lead 
effects are summarized by the many-spin Kondo Hamil- 
tonian. One effect is to enhance the ferromagnetic or- 
der. This result is important to make spintronic circuits 
by connecting leads in nanodevices. It is also promi- 
nent that a new peak appears in the specific heat but 
not in the susceptibility for small N. The peak position 
is Tk = (Jk/2J)Tc for the zigzag graphene nanoribbon 
lead. The energy is found to decrease around the peak 
position, and the entropy is lowered by factor /cb log 2 in 
the zero-temperature limit, indicating the Kondo effect 
due to a Kondo interaction between electrons in the lead 
and the nanodisk. 

We have proposed some applications of nanodisks to 
nanodevices. Being a ferromagnet, it can be used as a 
spin filter. Namely, only electrons with spin parallel to 
the spin of the nanodisk can pass through it. Addition- 
ally, it has a novel feature that it is not a rigid ferromag- 
net. The incoming spin-polarized current can rotate the 
nanodisk spin itself. Combining the advantages of both 
these properties, we have proposed a rich variety of spin- 
tronic devices, such as spin memory, spin amplifier, spin 
diode, spin valve and spin-field-effect transistor. Further- 
more, we have proposed some spin logic gates such as spin 
XNOR gate and XOR gate. Graphene nanodisks could 
well be basic components of future nanoelectronic and 
spintronic devices. 
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